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Quantum effect of one-dimensional photonic crystal 
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In this paper, we have studied the quantum transmission characteristics of one-dimensional pho¬ 
tonic crystal with and without defect layer by the quantum theory approach, and compared the 
calculation results of classical with quantum theory. We have found some new quantum effects in 
the one-dimensional photonic crystal. When the incident angle 9 = 0, there is no quantum effect. 
When the incident angle 9 ^ 0, we find there are obvious quantum effect with the incident angle 
increase. At the incident angle 9^0, there are also quantum effect with the change of thickness 
and refractive indexes of medium B or A. We think the new quantum effect of photonic crystal 
shall help us to design optical devices. 
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1. Introduction 

In 1987, E. Yablonovitch and S. John had pointed out that the behavior of photons. It can be changed when 
propagating in the material with periodical dielectric constant, and termed such material Photonic Crystal 
[1, 2]. Photonic crystal important characteristics are: Photon Band Gap, defect states, Light Localization 
and so on. These characteristics make it able to control photons, so it may be used to manufacture some 
high performance devices which have completely new principles or can not be manufactured before, such as 
high-efficiency semiconductor lasers, right emitting diodes, wave guides, optical filters, high-Q resonators, 
antennas, frequency-selective surface, optical wave guides and sharp bends [3, 4], WDM-devices [5, 6], 
splitters and combiners [7]. optical limiters and amplifiers [8-10]. The research on photonic crystals will 
promote its application and development on integrated photoelectron devices and optical communication. 
To investigate the structure and characteristics of band gap, there are many methods to analyze Photonic 
crystals including the plane-wave expansion method [11], Greens function method, finite-difference time- 
domain method [12-14] and transfer matrix method [15-20]. All of methods are come from classical Maxwell 
equations. In Refs. [21, 22], we have firstly studied the the quantum transmission characteristics of one- 
dimensional photonic crystal by the quantum theory approach, in which we have only considered the incident 
angle is zero, i.e., vertical incidence, we have found the classical and quantum transmission characteristics are 
the completely same, i.e., there is not quantum effect in one-dimensional photonic crystal. In this paper, we 
have studied the quantum transmission characteristics of one-dimensional photonic crystal when the incident 
angle is an arbitrary angle, i.e., non-vertical incidence, we find there are obvious quantum effect with the 
incident angle increase. At the incident angle 0^0, there are also obvious quantum effect with the change 
of thickness and refractive indexes of medium B or A. Otherwise, we have considered the effect of defect 
layer on the quantum transmission characteristics. When the incident angle 9 = 0, there is also not quantum 
effect. When the incident angle 9 ^ 0, with the incident angle increase, there are obvious quantum effect for 
the one-dimensional photonic crystal with defect layer. 

2. The quantum transmissivity 

In Refs. [23, 24], with the quantum theory approach, we have studied one-dimensional photonic crystal 
quantum transmissison characteristic when the incident photon is vertical incidence, i.e., the incident angle 
9 = 0. In the paper, we shall study one-dinrensional photonic crystal quantum transmissison characteristic 
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FIG. 1: The quantum structure of one-dimensional photonic crystal. 


when the incident photon is non-vertical incidence, i.e., the incident angle 9^0. 
The quantum wave equations of photon in medium is [23, 24] 


Vxf = 


E-V 

ch 





(1) 


where w is angle frequency of photon, c is the velocity of photon, ijj is the wave function of photon and n is 
refractive index of medium. 

The incident light, reflected light and transmission light are in the xz plane, the incident angle is 9 , which 
are shown in Fig. 1. The wave vectors A'/ and Kr of incident and reflection photon are 

Ki = Kj + K z k , I< R = -Kj + K z k (2) 


where K x = AT cos 0, K z = A'sin# are the wave vectors in x and z direction, and AT = ^, the w is the 
incident photon angle frequency. The wave functions ijii and i/jr of incident and reflection photon can be 
written as 


■ifj = F x e i( - K * x+K * z) i + F v e i( - K * x+K * z) j + F z e i{K * x+K * z) k, 
tf R = f^~ k * x+k * z) i + Fye i{ - K * x+K * z )j + F' z e*- K * x+K * z) k, 


( 3 ) 

( 4 ) 


where F x (F x ), F y (F y ) and F z (F z ) are the amplitudes of incident (reflected) photon wave functions in the 
x , y and z directions. 

Substituting ipi into Eq. (1), we have 


i 


k 


_d_ _d_ d_ 

dx dy dz 

jj^ e i(K x x+K z z) jp^ e i(K x x+K z z) jp^ e i(K x x-\-K z z) 


= -(F x e i( ~ K * x+K * z) i + F y e i{K * x+K * z) j + F z e i{K * x+K * z) k), 


dr 1 pi(K x x-\-K Z z) d jp* ( - ) i{K x x-\-K z z') uj jp „i(K x x-\-K Z z) 

dy rze dz r V t ' ~ c ’ 


JL F p i(K x x+K z z) I _d_jp p i(K x x+K z z) _ u jp p i(K x x+K z z) 
dx rze ^ dz rxC ~ c r V e ’ 


d jp pi^KxX+Kzz) _ _d_ jp ~i(K x x+K z z) _ pi{K x x-\-K z z) 

dx r y e dy rxt " — c rze ’ 


( 5 ) 


( 6 ) 
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From Eq. (7), we get 


-iK z F v = ±F X = KF X , 

- iK x F z + iK z F x = ^F y = KF y , 


iK x F v = ±F Z = KF Z , 


F -JJSlp 
tx — K ty, 


- iK x F z = (K -~^)Fy, 


K 2 x + K'l = K 2 . 


The incident current is [23, 24] 


Jj = icnjjj x 

r K 2 -I\ 2 ^ K z , l9 

= ^-KKf i+ -£ k]WA - (U) 

We find the incident current J is related to the amplitude F y of y component. So, we should only consider 
the j component wave function of photon in the following calculation. 

Firstly, we study the transfer matrices in the first period. The wave function of photon in medium A is 

= [Ale^ KAxX+KA ‘ z) + A2e i{ - KA:cX+KA ‘ z) }i 

-{-[A 1 ^ K AxX + K A zZ) + A , 1 e i (-KAa;X+KAzZ)l J 

■j-^A 1 e i ( KAxX+KAzZ ^ + A ,1 e i ^~ KAxX+KAzZ ^]k ( 12 ) 

where A x (A' x ), A l y ( Ay ) and Al {A' z ) are the amplitudes of incident (reflected) photon wave functions in 
the x, y and z directions. 

Substituting into Eq. (1), we have 

- ^-{A 1 y e i{KAxX+KA ^ z) + A£e*- KAmX+KA **)) = ^n A {Ale i(I<AxX+KA * z) + A^e i{ ~ KAxX+KA * z) ), ( 13 ) 

r) 8 

_ —_(A 1 qAKaxX+KazZ) _)_ A'^e^— K A xX+KazZ)\ _|_ _lj^^i(K Ax X+KAzZ) _|_ K A xX+K Az z)\ 

8x y z z dz x x 

= -n A [A 1 e i{KA * x+KA * z '> + A'y e i< '- KAxX+KAzZ ' ) ], ( 14 ) 


O (A l c i(K A * X +K Az z) _|_ ^/1 ^(-KaxX + KazZ)^ _ e i(KAxX + KAzZ) _|_ A ' 1 e^~ K AxX + Ka Z z) j 


From Eq. (13) to (15), we obtain 


Al _ _ Al 

Al 1 At 1 At 1 ' 

■™-y -™-x -^z 


( 16 ) 



In the incidence zone, the total wave function is 
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iptot(x,y,z) = i/> I (x,y,z)+ip R (x,y,z) 

_ (F x e i( ' KxX+KzZ> + F' e i ^ KxX+KzZ> )i 
+ (F y e i ( KxX+KzZ ) + f 1 e i ^~ KxX+KzZ ^)j 

+ ( Fze i( KxX + K * z ) + F' z e i{ - KxX+K * z) )k, (17) 

In the following, we should use the condition of wave function and its derivative continuation at interface of 
two mediums. 

(1) At x = 0, by the continuation of wave functions tytotix, y , z) and ip A (x, y, z), we have 

(F x + F' x )e iKzZ i + {F y + F' y )e iKzZ j + {F z + F' z )e iKzZ k 

= (A 1 , 1 + A' x 1 )e iKzZ i + (Ay + Ay )e iKzZ j + (A\ + A?)e iK ”k, (18) 

(2) At x = 0, by the derivative continuation of wave functions 4>tot(x,y, z) and ip A (x,y, z), we have 

%K X (F X - F' x )e iKzZ i + iK x (F v - F^)e iKzZ j 
+iK x (F z - F’ z )e iKzZ k 

= iK Ax (A x - A2)e iKAzZ i + iK Ax (Al - A^)e lKAzZ ] 

+iK Ax (A\-A*)e iKAzZ k, (19) 

with Eqs. (18) and (19), we get the j component relations 

(F v + F')e iKzZ = (A 1 + Ay)e iKAzZ , (20) 


K x (F y - F')e iKzZ = K Ax (A\ - Ay)e iKA * z , 


by Eqs. (20) and (21), we obtain 


K z = K Az , 


A-l + Ay — Fy + Fy 

A 1 - A ' 1 — Kx (F — F 1 ) ’ 
A y A v ~ k Ax \ r y r y> 


by Eq. (23), we can obtain the matrix form of A 1 A' 1 and F y , F' 


2 1 - 


= M\ 


L y \ 

F' ) ’ 


1 _ 

Ax 


where 


Mi = 


1 / 1 + 


2 1 - 


(25) 



is the transfer matrix of medium A in the first period, and the Kax is 
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K Ax = \j K‘a - K\ z = yjK\ - K 2 = sjK\ - K 2 sin 2 6, (26) 

where Ka = —n a , n a is the refractive indexes of medium A, and the transfer matrix M\ can be written as 


M\ = 


2 1 - 


The wave function of photon in medium B is 

ft B (x,y,z) = (Ble* KB * x+KB ‘ x) + B'^e i{ ~ KBxX+KBzZ) )i 

_|_ i(KsxX+KBzz) _|_ KbxX+Kb z z)^j 

+ (Ble KKBxX+KB * z) + B' z 1 e i{ ~ KBxX+KBzZ) )k, (28) 

where B\ ( B'£), By (By) and Bl (B' z ) are the amplitudes of incident (reflected) photon wave functions in 
the x, y and z directions. 

(3) At x = a, by the continuation of j component wave functions ^}\(x , y, z), ip B (x, y, z) and their deriva¬ 
tive, we have 

(A 1 + A'y)e iKAzZ + A\e iKAxa + A*e~ iKA * a = (By + B' 1 )e iKBzZ + Ble iI<Bxa + B''e~ iKB * a , (29) 


iI<Ax(A 1 e lKAxa+iKAzZ - A'}e- iKA * a+iKA * z ) = iK Bx (B 1 e 


il p iKBxa+iKBzZ _ D'l — iKB x a+iK B zZ 


with Eq. (29), we have 


(A 1 + A'he iKAzZ - (By + B'he iKBzZ = 0, 


with Eq. (30), we get 


A\e iKAza + A!}e~ iKA * a - (Ble iKBxa + B' y l e ~ iKBxa ) = 0, 


K A x(AI - Ay)e iIiAzZ - K Bx (B\ - B'Ae iKBzZ = 0, 


), (30) 


K Ax (A. 


l„iK Ax a \i\ 0 —iK Ax a 


- A'^e- zKAxa ) - K Bx (B. 


lJK Bx a + B 'l e -iK B xa^ = ^ 


with Eqs. (31) and (33), we obtain 


The Eqs. (32) and (34) can be written as 


I<Az = K Bz 


{ ^l e iK Ax a _|_ j^l e ~iK Ax a _ B 1 eiK Bx a , B ' 1 e~ iK Bxa 

K Ax (Ale iKA *“ - A£e~ iKAx “) = K Ex (Ble iKBxa - B£e~ iKB * a ) ‘ 

By Eq. (36), we can obtain the matrix form of By, By and Ay, Ay 

( B l\ _ 1 ( (1 + K Ax /K Bx )e i ^ KAx ~ KBx ' )a (1 -KAx/K Bx )e-^ KBx+KAx ^ \ ( A\ 


2 V (1 - K Ax /K Bx )e^ KAx + KB ^ a (1 + K Ax /K Bx )e i ( KBx - KAx > 
Mi (| ) , 



where 
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1 /( 1 + K Ax /K Bx )e^ K ^- K ^ a (1 - K Ax /K Bx )e-^ K ^+ KA ^ a 


1\/T 1 _ 

R " 2 V (1 - K Ax /K Bx )e i ( KA *+ K ^ a (1 + K Ax /K Bx )e^ K ^- K ^ a 

is the transfer matrix of medium B in the first period, and the K Bx is 


( 38 ) 


K Bx = JK%-K* sin 2 


(39) 


where K B = —rib, rib is the refractive indexes of medium B. 

Secondly, we use the similar approach can obtain the transfer matrices M A and M| of media A and B in 
the second period, they are 


M 2 _ 1 ( (1 + K B /K A )e i{KB ~ KA){a+b) (1 - K B /K A )e-^ KA+KB){a+b) 


and 


2\{1- K B /K A )e i ^ KA+KB ^ a+b '> (1 + K B /K A )e< KA - KB ^ a+b ^ 


2 _ 1 ( (1 + K A /K B )e i( - KA ~ KB ^ 2a+b '> (1 - K A /K B )e-^ KA+KB ^ 2a+b ^ 

Mb ~ 2 V (1 - K A /K B )e i( ~ KA+KB ^ 2a+b '> (1 + K a /K b ) e *(^B-^)( 2 “+ fe ) 


(40) 


(41) 


Finally, we can give the transfer matrices and Mg of media A and B in the N-th period, they are 


Ma =; 


1 ( (1 + K B /K A )e i{ - KB ~ KA ^ N - 1 ^ a+b ^ (1 - K B /K A )e - i< ~ KA+KB '> ( ~ N - 1 '> ( ~ a+b > 


2 V (1 - K B /K A )e i ^ KA + KB ^ N - 1 ^ a + b '> (1 + K B /K A )e i ^ KA ~ KB ^ N - 1 ^ a+b '> 


and 


Mr = 5 


1 / (1 + K A /K B )e i( ~ KA - KB '>( N ( a+b ')-V (1 - K A / _ft' B ) e -i(^+iCB)(JV(a+6)-6) 


2 V (1 — K a /K B ) e^ KA+KB ^ N ^ a+b) - b) (1 + K A /K B )e i(KB ~ KA)(N(a+b) - b) 

With the transform matrices, we can give their relations: 

(a) The representation of the first period transform matrices are 


A l \ _ M i ( Fv 
<) ~ A \F' y 


(42) 


(43) 


(44) 


b'Z)- M b \K 


= MhM\ (f V ')= Ml ( p] 


(45) 


(b) The representation of the second period transform matrices are 


A l\-M 2 ( B y 
<-)- A \B* 


- ( Fy i)= ( F p y, 


p2 

J2 I = M 
B v 


b ( §2 ) = M 2 B M 2 A MlM\ ( ^ ) = M 2 M l 


(46) 

(47) 


(c) Similarly, the representation of the N-th period transform matrices are 


A n 

S±y 

A' N 
A y 


= M A M b ~ 1 M a ~ 1 ■ ■ ■ M 2 A M l B M\ [ p'/J = m a m n_1 ■ ■ ■ m 2 m 1 ^ p v , 


(48) 



tjN 

B y | — 
R IN ~ 

B V 


-Uj'; A/^A/,'; 1 A/ “' • • ■ j = hl N M N ~ l ■ ■ • M 2 M l ( (j j = A/ 


$ I (®) 


where 


M = MM 


N A fN-l 


m 2 m 1 = f TOl W2 

m3 WI4 


(50) 


is the total transform matrix of N period, and M 1 = MgM\ is the first period transform matrix, M 2 = 
M%M\ is the second period transform matrix, and M N = Mg M% is the N-th period transform matrix. 
The wave function of N-th period in medium B is 


$B{x,y,z) = (B* e i(KB * x+KB * z) + B X N e i{ - KB * x+KB * z) ji 
+ (B n e i ^ KBxX+KBzZ ' ) + b n e i ^~ KBxX+KBzZ ^)j 

+ (B N (^(KbxX+Kbzz) +B 'N e i(-K Bx x+K Bz z)^ 

In FIG. 1, the transmission wave function is 

$ D (x, y, z) = D x e i ^ x+Kzz) i + D v ei( K * x + K * z )j + D z e i(K * x+K * z) k. 


(51) 


(52) 


(4) At x = N(a + b), by the j component continuation of wave functions ipg(x, y, z) and i/>d(x, y, z), we have 


B N e i(K Bx N(a+b)+K Bz z) + B 'N e i(-K Bx N(a+b)+K Bz z) _ ]j^ e i(K x N(a+b)+K z z) ^ 

since the Eq. (53) is an equation for an arbitrary variable z , we have 

B N e iK Bx N(a+b) _|_ B 'N e -iK Bx N(a+b) _ jj e iK x N(a+b) 

with Eqs. (49) and (50), the Eq. (54) can be written as 

(mi F v + rri 2 Fy)e iKBxN< ' a+b ^ + ( m 3 F y + m 4 F^e- iKB * N ( a+b) = D y e iK * N ( a+b \ 


(53) 


(54) 


(55) 


(5) At x = N(a+b ), by the j component derivative continuation of wave functions ipg (x, y, z) and i[>d(x , y, z), 
we have 


K Bx B”e iK *- N W> - K Bx B y N e- iKB * N ( a+ V = K x D y e iK ^ a+b \ 
with Eqs. (49) and (50), the Eq. (56) can be written as 


(56) 


^(mi F y + m 2 F')e iKB * N ( a+b) - ^( m 3 F y + rruF')e- iKB * N ( a+ V = D v e lK ^ a+b \ (57) 

K r y K r y 


By Eqs. (55) and (57), we can obtain 


F' 

v 


mi{K x - K Bx )e iKB * N ^ +b ') + m 3 (K x + K Bx )e~ iKB - N ^ 
m 2 {KBx — K x )e lKBxN ( a+b '> — m^(K x + Kbx)^~ iKBxN ( a+b ^ 


(58) 


t = 


£y 

F, 


= (mj + + (m 3 + m ^)e~ i ^+ K ^ a + b \ 


F, 


and the quantum transmissivity T is 


T=\tf 


(59) 


(60) 





(b) 0=0 classical ®/© 0 


FIG. 2: The quantum and classical transmissivity of incident angle 9 = 0 . (a) quantum transmissivity (b) classical 
transmissivity. 


5. Numerical result 

In this section, we report our numerical results of quantum transmissivity. The main parameters are: For 
the medium A, its refractive indexes is n a = 1.45, and thickness is a = 267 nm. For the medium B , its 
refractive indexes is rib = 3.59, and thickness is b = 108 nm.. The central frequency is u o = 1.216 x 10 15 Hz, 
and the period number is N = 16. With Eqs. (58)-(60), we can calculate the quantum transmissivity. 
In FIG. 2, we calculate the quantum and classical transmissivity when the incident angle 9 = 0 , FIG. 2 
(a) and (b) are quantum and classical transmissivity, respectively. Comparing FIG. 2 (a) and (b), we find 
the quantum transmissivity and classical transmissivity are the completely same, i.e., when the incident 
angle 9 = 0, there is not quantum effect in one-dimensional photonic crystal. In FIG. 3, we calculate the 
quantum and classical transmissivity when the incident angle 9 are respectively. From FIG. 

3 (a) to (d), they are classical transmissivity. From FIG. 3 (e) to (h), they are quantum transmissivity. 
We can obtain the following results: (1) With the incident angle increase, the forbidden bands width are 
unchanged and positions red shift for the classical transmissivity. (2) With the incident angle increase, the 
forbidden bands become widened and positions blue shift for the quantum transmissivity. (3) For the same 
incident angle 9, the quantum forbidden bands are wider than the classical forbidden bands, and the quantum 
forbidden bands positions blue shift. (4) When the incident angle increase, the quantum effect become more 
remarkable. In FIG. 4, at the incident angle 9 = we calculate the quantum and classical transmissivity 
when the thickness b of medium B are 108nm and 158mn, respectively. FIG. 4 (a) and (b) are the classical 
transmissivity, and FIG. 4 (c) and (d) are the quantum transmissivity. We can obtain the following results: 
(1) With the thickness b increase, the forbidden bands numbers increase and positions red shift for the 
classical and quantum transmissivity. (2) For the same thickness b, the quantum forbidden bands are wider 
than the classical forbidden bands, and the quantum forbidden bands positions blue shift relative to the 
classical forbidden bands. In FIG. 5, at the incident angle 9 = we calculate the quantum and classical 
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FIG. 3: The quantum and classical transmissivity for different incident angle 9 . (a)-(d) classical transmissivity, (e)-(h) 
quantum transmissivity. 


transmissivity when the refractive indexes rib of medium B are 3.59 and 4.09, respectively. FIG. 5 (a) and 
(b) are the classical transmissivity, and FIG. 5 (c) and (d) are the quantum transmissivity. We can obtain 
the following results: (1) With the refractive indexes rib increase, the forbidden positions red shift for the 
classical and quantum transmissivity. (2) For the same refractive indexes rib, the quantum forbidden bands 
are wider than the classical forbidden bands, and the quantum forbidden bands positions blue shift relative 
to the classical forbidden bands. In FIG. 6, we calculate the quantum and classical transmissivity with defect 
layer, and the incident angle 9 = 0, FIG. 6 (a) and (b) are quantum and classical transmissivity, respectively. 
Comparing FIG. 6 (a) and (b), we find the classical and quantum transmissivity are identical, i.e., when 
the incident angle 9 = 0, there is not quantum effect in one-dimensional photonic crystal with defect layer. 
In FIG. 7, we calculate the classical and quantum transmissivity with defect layer, and the incident angle 
are |, f, respectively, and the structure is ( AB) 8 D(AB ) 8 . From FIG. 7 (a) to (c), they are classical 
transmissivity. From FIG. 7 (d) to (f), they are quantum transmissivity. We can obtain the following results: 
(1) With the incident angle increase, the forbidden bands width and defect model intensity are unchanged, 
but the forbidden bands and defect model positions red shift for the classical transmissivity. (2) With the 







10 



FIG. 4: The quantum and classical transmissivity of incident angle 9 = for different thickness b. (a)-(b) classical 
transmissivity, (c)-(d) quantum transmissivity. 


incident angle increase, the quantum forbidden bands become widened, the defect model intensity weaken, 
and the forbidden bands and defect model positions blue shift relative to the classical transmissivity. (3) For 
the same incident angle, the quantum forbidden bands are wider and defect model intensity weaker than the 
classical forbidden bands, the defect model and the quantum forbidden bands positions blue shift relative to 
the classical. (4) When the incident angle increase, the quantum effect become more remarkable. 
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FIG. 5: The quantum and classical transmissivity of incident angle 0 = for different refractive indexes rib■ (a)-(b) 
classical transmissivity, (c)-(d) quantum transmissivity. 


6. Conclusion 

In summary, we have studied the quantum transmission characteristics of one-dimensional photonic crystal 
by the quantum theory approach, and compared the calculation results of classical with quantum theory. 
We have found some quantum effects in one-dimensional photonic crystal. When the incident angle 0 = 0, 
i.e., vertical incidence, the classical and quantum transmission characteristics are the completely same, i.e., 
there is not quantum effect in one-dinrensional photonic crystal. When the incident angle 9 ^ 0, we find 
there are obvious quantum effect with the incident angle increase. At the incident angle 9 ^ 0, there are also 
obvious quantum effect with the change of thickness and refractive indexes of medium B or A. Otherwise, 
we have considered the effect of defect layer on the quantum transmission characteristics. When the incident 
angle 0 = 0, there is also not quantum effect, and when the incident angle 0^0, with the incident angle 
increase, there are obvious quantum effect for the one-dimensional photonic crystal with defect layer. The 
new quantum effect of photonic crystal shall help us to design optical devices. 
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FIG. 7: The quantum and classical transmissivity for different incident angle 9 with defect layer, the structure is 
(AB) 8 D(AB) 8 . (a)-(c) classical transmissivity, (d)-(f) quantum transmissivity. 








